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A6sfracf-Bracewell has proposed the Discrete Hartley Transform (DHT) as a substi- 
tute for the Discrete Fourier Transform (DFT), particularly as a means of convolution 
(R. N. Bracewell, "Discrete Hartley transform,” Journal of the Optical Society of Amer- 
ica, vol. 73, pp. 1832-1835, 1983 ). Here we show that the most natural extension of the 
DHT to two dimensions fails to be separable in the two dimensions, and is therefore 
inefficient. We consider an alternative separable form, and derive a corresponding con- 
volution theorem. We also argue that the DHT is unlikely to provide faster convolution 
than the DFT. 
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1. Introduction 


In a series of recent papers Bracewell has drawn attention to the Hartley transform 
as a substitute for the much more widely known Fourier transform. [1, 2, 3] He has 
introduced a discrete Hartley transform (DHT) as an analog to the discrete Fourier 
transform (DFT), and has developed a fast algorithm (FHT) to compete with the Fast 
Fourier Transform (FFT). The purpose of this paper is to point out certain problems 
encountered in the extension of the DHT from one dimension (ID) to two dimensions 
(2D). There are two possible ways in which this extension can be made. Each has its 
virtues, but only one is practical. This choice requires a new convolution theorem, 
which we provide. However, even with this theorem, the 2D DHT appears not to pro- 
vide a computational advantage over the DFT. 


2. The Discrete Fourier Transform (DFT) 

Consider a discrete sequence p (ar ), x = 0, • • • ,N -1. The forward DFT of p (x ) 
is given by 


N - 1 

P F (« ) = E P i x ) ex P(“* Znux /N ) (1) 

x =0 

and the inverse DFT, 


P(*) = 4; E Pf ( u ) exp(t 27m; /N ) (2) 

•* v u =0 

In our notation we depart slightly from Bracewell, who places the scalar l/N in the for- 
ward transform in both the DFT and the DHT. This differs from general practice, 
[4, 5, 6] and is less efficient if forward transforms are more common than reverse 
transforms. We therefore place the scalar in the inverse transform in our definitions of 
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the DFT, DHT, and SDHT (defined below). 

3. Symmetry in the DFT of a real sequence 

Although the DHT is defined for both real and complex sequences, its practical 
value arises from the way in which it takes advantage of the symmetry in the DFT of a 
real sequence. Therefore in what follows we consider only real input sequences. A real 
sequence, defined by N real numbers, has a DFT that contains N complex entries, or 
2 N real numbers. However, only N of the 2 N real transform values are independent 
because the transform is Hermitian, that is, its real part is even and its imaginary part 
is odd. The Hartley transform is one way of taking advantage of this redundancy. 
Note that the sum (or difference) of an odd and an even sequence can be easily broken 
apart again into its odd and even parts. Thus because real and imaginary parts of the 
DFT are even and odd respectively, we can add them together with the assurance that 
they can be broken apart to recover the full complex DFT when desired. We therefore 
define the ID DHT of a real sequence p (x ) as the real part minus the imaginary part of 
the DFT, 

Ph (« ) = RePp (« ) - ImP/r (u ) (3) 

Expanding the right side by means of Eq. 1, 

N - 1 

Ph{ u )= E P (* ) I cos2 nut /N + sin2^uar /N ] (4) 

s =0 

Ph (« ) = E P ( x ) cas(2^«2 /N ) (5) 

x =0 

The sequence cas defined by 

cas(2? tux /N ) = cos (2t tux /N ) + sin(2^uar /N ) (6) 

therefore serves as the basis sequence for the DHT. Note that the DHT may be 
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expressed as the sum of its even and odd parts 

P H ( w ) = P HE ( u ) + P HO ( u ) ( 7 ) 

where 

P He( u ) = *|[ P h( u )+ P h(~ u ) 1 P HO ( M ) = *1 1 P H ( u )~ P H (-« ) ] (8) 

Throughout the text we interpret all sequence arguments modulo the length of the 

sequence, so for example if P (u ) is a sequence of length N then P (~l)=/ > (N - 1). 

Note that the DFT is easily reconstituted from the even and odd parts of the DHT 

P F («) = Phe (« ) - <Pho (« ) • ( 9 ) 

The inverse transform is 

i N - 1 

p (*) = “Ar S P h( u ) c^{2t:ux/N) (10) 

iV « =o 

The principal advantage of the DHT over the DFT, as pointed out by Bracewell, is that 
it transforms a real input into a real output. Furthermore, the convolution theorem can 
be expressed in terms of the DHT. While this expression involves four sequence pro- 
ducts, rather than the single product of the DFT convolution theorem, the former are 
products of real quantities and will therefore require on the order of 1/4 the number of 
real operations. Furthermore, if either of the inputs to the convolution is even, as is 
quite common in filtering applications, the result reduces to a single product of real 
sequences, requiring on the order of 1/4 as many operations as convolution via the DFT. 

4 . The DHT in two dimensions 

One of the most frequent applications of the DFT is in image processing, where it is 
frequently used to implement convolution. For the DHT to be of use in image process- 
ing it must be extensible to two (and more) dimensions, and these extensions must carry 


4 



with them a simple form of the convolution theorem. In one of his papers, Bracewell 
has offered a definition of the 2D DHT. [2] But this definition does not agree with the 
general ID definition he proposes as embodied in Eq. (9), or with the convolution 
theorem he has offered (Ref. 2, p. 1834). 

There are two natural candidates for the extension of the DHT to two dimensions. 
The basis functions of these two possibilities are 

cas[ 2 n(tix /N +vy /M ) ] and cas(2;mr /N )cas(2 xvy /M ) . (11) 

In what follows, we will consider the advantages of each alternative. 

First, which definition is consistent with the definition of the Hartley as the 
difference of even and odd parts of the DFT? Recall the observation that motivated the 
construction of the ID DHT, namely that the DFT of a real input is Hermitian. Let 
p (x ,y ) be a real 2D sequence (image) with x = 0,...,iV-l and y = 0,...,A/-1. In two 
dimensions, odd-ness and even-ness are defined with respect to reflections about both 
axes: 

Even {p (x ,y )} = p E (x ,y ) = | [ p (x ,y) + p (-x -y ) ] (12) 

Odd {p (x ,y )} = vo (* ,y ) = | [ p (* ,y ) - p (-* -y ) ] • (13) 

Under these definitions, the real part of the DFT of a real input is even, and the ima- 
ginary part is odd. We can again define the DHT to be the even part minus the odd 
part of the DFT, 

P H («,«)= R eP F (« ,v ) - Im P F (« ,v ) (14) 

= Phe ( tt ) + Pho ( u > v ) • f 1 ^) 

In this case the DFT, as for ID, is given by 
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P F (u ,v ) = P HE (u ,v ) - iP„o (ti ,v ) . (16) 

To discover the basis sequences of the 2D DHT we expand Eq. (15). Let 

a = 2nux /N and 6 = 2wvy /M . Then 

N - 1 M - 1 

Ph (« ) = E E P (* )[ cos ( a + b ) + sin(a +6 ) ] (17) 

x =0 y =0 

N-l M - 1 

= E E p(*,y)cas(a+6) . (18) 

x =0 y =0 

Thus we see that a definition of the 2D DHT that preserves the symmetries of the ID 
transform leads to the first of the two possible basis functions introduced above. 

This definition, while it forms a natural extension of the ID transform, and leads to 
the convolution theorem offered by Bracewell, [2] has a serious drawback. Unlike the 2D 

’ i 

DFT, it is not separable in « and v . Fast algorithms for the 2D DFT rely upon separ- 
ability, and a straightforward extension of the FHT to 2D would also require separabil- 
ity. 

* 

5. The Separable 2D Hartley Transform 

As a solution to this dilemma, we consider an alternative definition of the 2D DHT, 
which we will call the Separable Discrete Hartley Transform (SDHT), defined by 

N-l M -1 

Psi u > v )= E E P (x ,y ) cas(2ff«x /N ) cas(2xvy /M ) (19) 

x =0 y =0 

and the inverse transform 

j N-l M-l 

P( x ’?) = UD E E P S {^ iv) cas(2irux /N)cas(27rvy /M) . (20) 

ism II =0 t =0 

This is essentially the 2D DHT as given by Bracewell (only the position of the scale fac- 
tor is different). [2] Note that both transforms are separable in x and y , and can be 
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executed by first applying a ID DHT to each row and then a ID DHT to each column 
(or vice versa). Since this transform no longer obeys Eq. (9), and does not lead to the 
convolution theorem derieved therefrom, we must discover the relationship between the 
SDHT and the DFT, and derive a convolution theorem for the SDHT. 

0. Relation between SDHT and DFT 

In the SDHT we have gained separability but have lost the simple relationships 
between the DHT and the DFT, namely, DHT = ReDFT - ImDFT and 
DFT = Even DHT - i Odd DHT . However, as we show below, a similarly simple 
relation exists for the SDHT and the DFT. We begin by expanding the cas terms in the 
definition of the SDHT, 

N- 1 M- 1 

P s (« ,v ) = 2 J] p (x ,y ) (cosfl 4- sina ) (cos6 4- sin6 ) (21) 

x =0 y =0 

N-i M-l 

— Yj p (x ,y )(cos a cos6 4- sina sinft + cosa sin6 + sina cosft ) (22) 

x =0 y =0 

N- 1 M-l 

= E E p (* »y ) [ cos ( a ~ b ) + sin ( a + 6 ) 1 (23) 

* =0 y =0 
N- 1 M-l 

= E E p(x ,y){Re[cos(a-b)-ism(a-b)\ (24) 

x =0 y =0 

- Im[ cos(a +b )-i sin(a 4-6 ) ] } 

= R eP F (u -v ) - lmP F (u ,v ) . (25) 

Thus the SDHT is equal to the real part of the DFT reflected about the u axis, minus 

the imaginary part of the DFT. As might be imagined, it does not matter whether the 

reflection is about the ti or# axis. 
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( 26 ) 


Going in the other direction, we can expand the expression for the DFT, 

Pp (« ,v ) = RePp (« ,v ) + i lmP F (u ,v ) 

= — { [ RePp (« ,v ) -f RePp (-« ,-v ) - Im P F (u ,-v ) - Im P F (-« ,v ) ] (27) 

- » [ R ePp (-« ,v ) - RePp (u ,-v ) + ImPp (v ,v ) - ImPp (-u ,-v ) ] } 

= 1 p s («.-») + Ps (-«.») 1 - Ij, I Ps («.») + Ps (-« -w )Md 1 (28) 

= PsE (« ) ~ fPsO {«,'>) ■ (29) 

7. Summary of relationships 

For future reference we summarize the relationships among the DFT, DHT, and 


SDHT: 

P s (« ,t> ) = P ho («,”) + P HE (« ,-v ) (30) 

P s (u ,v ) = RePp (u ,-v ) - ImPp (u ,v ) (31) 

P H (« ,v ) = P so (« ,v ) + P S E (« ,-v ) (32) 

P H (u ,v ) = RePp (u ,v ) - ImPp (u ,v ) (33) 

Pp (u ,v ) = P SE (u ,-v ) - iP so (u ,v ) (34) 

Pp (« ,v ) = P HE (« ,v ) - iP H0 (u ,v ) . (35) 

8. SDHT Convolution Theorem 


To be useful, the SDHT must yield a simple convolution theorem. We can derive 
this theorem from the preceding relations between DFT and SDHT. Let r (x ,y ) be the 
result of convolving the two real sequences p (x ,y ) and q (x ,y ), 

N - 1 M - 1 

r(z ) = p (* ,y ) * q (* ,y ) = £ £ p(*J) g • (36) 

«' =0 / =0 

From the convolution theorem for the DFT we have 
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( 37 ) 


R f («,») = P F (« ,# )Qf (u,v) . 

From Eq. (31) we can convert this to an expression for the SDHT, 

R s (« ,v ) = ReRp (« ,-v ) - lmR F (« ,v ) (38) 

= Re[ P F (« ,-v )Q F (ti ,-v ) ] - Im[ P F (« ,v )Q F (u ,v ) ] . (39) 

The next step is to substitute Eq. (34) for each DFT, 

= Re{ [ P SE (« ,v ) - i P so (* ,-v ) ][ Q se (« ,v ) - iQ S0 (u ,-v ) ] } (40) 

- Im{ [ P SE (« ,-v ) - iP S0 (« ,v ) ][ Q se (u ,-v ) - iQ S0 (« ,v ) ] } . 

Thus we arrive at a convolution theorem 

R s (« ,v ) = P SE (« ,v )Q 5£ (« ,v ) - P so (u ,-v )Q 50 (« ,-v ) (41) 

+ Pse (« ~ v ) Qso («»«) + Pso (« > v )Qse (« ~ v ) • 

9. Symmetrical inputs 

Bracewell notes that if one of the inputs to a convolution is even, then the DHT of 
the convolution reduces to the product of the DHTs. For the SDHT, examination of Eq. 
(41) shows that if p (x ,y ) is even, the SDHT of the convolution reduces to 

R s (« ,v ) = P SE (« ,v )Qse (« ) (42) 

+ Pso (« )Qso (« ~ v )+ p se (« ~ v )Qso (« .« ) • 

Thus only three products need to be computed. 

10. Complexity 

As noted above, the SDHT is one method of exploiting the symmetry of the DFT of 
a real sequence. Other well-known methods allow one to compute the DFT of a real 
sequence in essentially half the number of operations required for a complex sequence. 
(5] Let us call this sort of transform an RDFT. Furthermore, the result of an RDFT on 
a real sequence of length N is N /2 complex values. Following Bracewell, let us assume 


0 



for the moment that for a sequence of length N the number of operations required by 
SDHT and RDFT are equal. The convolution via the RDFT will require N /2 complex 
multiplies, or approximately 2 N real multiplies. Convolution via the SDHT will require 
4iV real multiplies in the general case. When one of the inputs is symmetrical, then as 
noted above only 3 N real multiplies are required. In neither case, however, does the 
SDHT confer an advantage over the DFT. 

11. Summary 

We have shown that the most natural extension of the DHT to two dimensions fails 
to be separable in the two dimensions. We have examined an alternative, separable 
form for the 2D DHT, originally proposed by Bracewell. We have determined its rela- 
tion to the DFT and have derived a corresponding convolution theorem. Finally, we 
note that the DHT is unlikely to provide faster convolution than the DHT in the general 
case. 
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